Pseudo-magnetic field distribution and pseudo-Landau 
levels in suspended graphene flakes 

M. Mucha-Kruczyhski*' a , V.I. FaPko a 

a Department of Physics, Lancaster University, Lancaster, LAI J^YB, UK 



Abstract 

Combining the tight-binding approximation and linear elasticity theory for a 
planar membrane, we investigate stretching of a graphene flake assuming that 
two opposite edges of the sample are clamped by the contacts. We show that, 
depending on the aspect ratio of the flake and its orientation, gapped states 
may form in the membrane in the vicinity of the contacts. This gap in the 
pre-contact region should be biggest for the armchair orientation of the flake 
and ^ < 1. 
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Over the past eight years since its successful isolation [1], graphene has 
attracted a lot of attention [2, 3, 4]. In particular, fabrication of high quality 
graphene samples suspended over the substrate and thus decoupled from its 
influence [5, 6] stimulated investigations of the closest vicinity of neutrality 
point in the electronic band structure looking for electron-electron interaction 
effects in both monolayer [7] and bilayer [8, 9, 10, 11] graphene. 

However, it is easy to imagine that such an atomically flat membrane placed 
over a ridge undergoes unintentional mechanical deformations during the fab- 
rication process or operation of the device. At the same time, it is well-known 
that for graphene systems lattice distortions are coupled to electronic degrees 
of freedom resulting in gauge fields resembling valley- dependent effective mag- 
netic field B acting on electrons in the deformed flake [12, 13, 14]. It has been 
shown that application of an electric field with a gate placed underneath the 
flake deforms it and leads to suppression of the conductance [15]. Similarly, in a 
four-contact configuration the interaction between the side contacts and charge 
carriers in graphene distorts the flake and can impede observation of quantum 
Hall plateaus [16]. 

In this paper, we combine the tight-binding approximation with the linear 
theory of elasticity and investigate distributions of the strain-induced pseudo- 
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magnetic field with relation to the flake orientation and its aspect ratio. We 
show that for the zig-zag orientation of the flake, a B = region connecting 
both contacts always exists in the middle of the flake. In contrast, for the 
armchair orientation, regions of nonzero pseudo-magnetic field are situated along 
the contacts and we highlight the possibility of strain-induced Landau level 
formation within those regions for samples with W < L. 

We model the suspended flake with a flat rectangular elastic membrane of 
length L and width W clamped at the edges by metallic contacts, as shown in 
Fig. 1(a). We denote the membrane coordinate system as xy and the crystalline 
coordinate system as x'y' . We choose the latter in such a way that for monolayer 
graphene the x' axis points along the armchair direction and call the angle 
between the x and x' axes. 

After the deformation, the total elongation of the membrane along x is dL 
(4p on each side), leading to an inhomogeneous distribution of strain. Here, 
we are interested in the effects caused by small strains (< 1%) and ignore 
any strain-induced wrinkling developing in the membrane. Then, the problem 
simplifies to that of stretching of a planar sheet. Within the linear elasticity 
theory [17, 18], the equations for displacement u(x,y) = (u x ,u y ) of a point 
(x, y) in the membrane take the form 

2>d xx u x + (1 - cr)dyyU x + (1 + (j)d X yU y = 0, 

2dyyUy + ( 1 ~ (l)d XX Uy + (l + (l)d X yU x = 0, 

with clamped boundary condition for the left and right edges and free boundary 
condition for the top and bottom edges, 

(u x (±±,y) = ±f i[<Td x u x +d y u y ] y=± m=0 
\u y (±%,y)=0 ' \[d x u y + d v u x ] y=±i v=0 ■ {> 

Above, a is the Poisson's ratio for graphene and in this work, we take a = 0.12 
[19, 20, 21]. 

However, microscopically the lattice is described in the crystalline coordi- 
nates x'y' . On this level and within the tight-binding picture, strain-induced 
changes in the distances between carbon atoms in the lattice lead to shifts of 
the energies of the electron 7r-orbitals, €2 P . This introduces the diagonal term 
^2 P ^ab y # u with e' 2p — q1^ b ~ 6eV/A[22] expressing the change of the on-site 
energy with growing bond length tab (for unperturbed lattice, tab — 1.46 A). 
At the same time, change in the bond lengths introduces asymmetry in the 
couplings for the three nearest neighbours surrounding any given atomic site. 
The form of the resulting addition to the electronic Hamiltonian written for the 
states in the vicinity of one of the Brillouin zone corners (valleys) is equivalent 
to the addition of a valley- dependent gauge vector potential [3], 
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The parameter 770 = d\nrAB ex P resses the change of the nearest neighbour 
coupling 70 with the change of the bond length tab- Its value 770 ~ — 3 can 
be estimated from the DFT calculations [22] and electron-phonon interaction 
constant measured by Raman spectroscopy [23]. The u x > and u y > are the com- 
ponents of the displacement in the crystalline coordinates, £ = ±1 identifies one 
of the two valleys and e is electron charge. A strong enough inhomogeneous 
strain may cause quantization of electronic states into Landau levels [24, 25, 26] 
and introduce gaps into the electronic spectrum. With the help of the angle 
6 expressing the connection between the membrane and crystalline coordinate 
systems, B can be written in the macroscopic coordinates xy as 
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B = o _ _ {sin3fl [{d xx - d yy ) u x -2d xy u y ] -cos30 [2d xy u x + (d xx -d yy )u y ]}. 



The problem of a planar membrane stretched like in Fig. 1(a) does not have 
a simple analytical solution and hence we applied the finite element method 
[27] to obtain displacements in membrane coordinates. Due to the need of the 
first and second derivatives of displacement, a nine-point rectangular element 
has been used. An example of the u x and u y components of the displacement 
for a membrane with length L = 1/im, aspect ratio ^ = 1 and longitudinal 
strain ^& = 1% is presented in Fig. 1(b) 1 . The first derivatives of the dis- 
placement which appear in the expression for the gauge potential *4o, Eq. (3), 
are associated with the components of the stress tensor. For planar problems 
treated within the linear theory of elasticity, stress singularities occur in sharp 
corners [28] due to the incompatibility of the boundary conditions at the cor- 
ner point and nonlinear theory is required to obtain physical results in their 
vicinities [29, 30]. Hence, within our model, both Aq and B diverge at the cor- 
ners of the flake. These divergencies are only important within several lattice 
constants away from the corner [31] and we ignore those regions in our consider- 
ations. Physically, the boundary conditions can be locally relaxed as the stress 
is limited by the treshold for sliding of the flake. 

Having obtained the first derivatives of the displacement, one can compute 

the strain-induced on-site scalar potential <&(x,y) — e2p ^ AjB V • u{x,y). How- 
ever, as opposed to the vector potential *4o, the scalar potential is screened 
by carriers in the flake. Detailed investigations of this screening indicate that 
linear response theory approximates its effects reasonably well [32] . Within this 
scheme, the induced carrier density 5n(x,y), is [32] 



5n(x, y) = T 



(5) 



where x(o) 18 the static zero-temperature Lindhard function of a homogeneous 
noninteracting massless Dirac fermion fluid, e(q) is the static random-phase- 



1 A mesh of 400 x 400 elements used in this example is a good representative of the mesh 
sizes used in this work. 



3 



approximation dielectric function and J r [. . .] denotes a Fourier transform. We 
have checked that for our problem Sn is negligible and thus the scalar potential 
can be ignored. 

In Fig. 2, pseudo-magnetic field B seen by electrons in the vicinity of the K+ 
valley (£ = 1) for a square membrane of length L = 2/im and ^ = 1% is shown 
in relation to the angle (the pseudo-magnetic field for electrons in the other 
valley has the opposite sign). As the angle is varied, the distribution of B changes 
drastically. In particular, for = [Fig. 2(a)], B is asymmetric about the x 
axis. It follows that electrons propagating along the y = line do not experience 
any perturbation associated with the pseudo-magnetic fields. However, as 
increases, this B — region in the middle of the flake is destroyed. The areas 
where the magnitude of B is significant shift from the centres of individual 
quarters of the flake towards the left and right edges [Fig. 2(b) and (c) for 
= and respectively] and for = J [Fig. 2(d)] B is concentrated along 
those edges and antisymmetric about the y axis. This configuration may severely 
impact the transport through the system if electronic states are quantized into 
Landau levels as a result of B. To inspect the latter possibility, we look at the 

magnitude of the 'magnetic length' A# associated with the field 6, Xb = \ \, 



as an indicator of the stability of the Landau level formation: if the pseudo- 
magnetic field B changes little at the distance comparable with local A#, the 
Landau level may form in that portion of the flake. For that purpose, three 
points on the flake, marked with black, green and orange dots, have been chosen 
in Fig. 2 and a circle with a radius corresponding to the magnetic length A# at 
those points have been drawn using respective colours, for each of the graphs 
in Fig. 2. We see that in the case of = |, quantization of electron states 
into Landau levels may be indeed expected along the left and right edges of 
the flake: the orange and black circles in Fig. 2(d) cover regions of reasonably 
smooth pseudo- magnetic field B ~ 0.1T. For monolayer graphene, that suggests 
a gap of the order of lOmeV between the zero-energy Landau level and the 
next quantized state. To compare, in bilayer graphene (within our model, the 
distribution of B is the same irrespectively of the number of graphene layers 
in the flake), the gap would be only of the order of 0.5meV. In contrast, no 
pronounced quantization of the electronic states should be expected for the 
flake orientation at = 0. 

We now investigate the dependence of the distribution of the pseudo-magnetic 
field B on the aspect ratio of the suspended flake. In Fig. 3, the strain-induced 
distribution of B is shown for flakes with aspect ratios ^ = 2 and ^ = 0.5. 
We consider two most characteristic orientations of the flake with respect to the 
x axis, the 'zig-zag' (6 = 0) and 'armchair' (9 = |) orientations. Again, we 
choose three points on the flake for which a circle of radius equal to the local 
magnetic length A# is drawn. For flakes in zig-zag orientation, a possibility of 
Landau level formation in some areas is indicated in the case of wide and short 
samples. However, due to symmetry, for such orientation the top and bottom 
regions of nonzero B are always separated by a region of B = 0. For the arm- 
chair orientation, the regions of significant B are located in the vicinity of the 
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contacts. Irrespectively of the flake orientation, for very wide or very narrow 
flakes, the centre of the flake is far enough from the corners to perceive the 
strain as homogeneous and such strain does not generate any pseudo-magnetic 
field [24, 25, 33]. 

In summary, we have investigated the distribution of the deformation-induced 
pseudo-magnetic field in suspended graphene flakes stretched along two oppo- 
site edges which remain clamped 2 . We show that details of the distribution of B 
depend on the microscopic orientation of the flake with respect to the contacts. 
In particular, the zig-zag orientation always possesses a B = region connecting 
both contacts in the middle of the flake. For the armchair orientation, similar 
region connects the free top and bottom edges of the flake and the pseudo- 
magnetic flake is concentrated in the vicinity of the contacts. We suggest that 
for the latter, Landau levels may form along the left and right edges due to the 
presence of £>, especially in the flakes with aspect ratios ^ < 1. Such situation 
would create within the flake regions with gapful electronic spectrum and we 
estimated the magnitude of this gap as ~ lOmeV for square monolayer graphene 
flakes with the side length of 2/im. 
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Figure 1: (a) Rectangular graphene membrane under tension. The red dashed and blue 
solid line depict the membrane before and after the deformation, respectively. In the inset, 
the relation between the membrane and crystalline coordinates xy and x'y' is shown, (b) 
Distribution of the u x and u y components of displacement for a square membrane. 
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Figure 2: Distribution of the pseudo-magnetic field B for electrons around the K+ valley as 
a function of the angle for a square membrane with L = 2/xm and ^ = 1%. Note that the 
distribution of B in the vicinity of the sharp corners of the flake is not shown and replaced 
with white squares. The black, green and orange circles show the magnetic length for three 
points in the membrane marked by a dot in the same colour. 
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Figure 3: Distribution of the pseudo-magnetic field B for electrons around the K+ valley for 
flakes with aspect ratios ^ = 2 (top) and ^ = 0.5 (bottom). Like previously, B in the 
vicinity of the corners is ignored. In all cases, ^ = 1%. Distributions for both 6 = and 
= ^ are shown. The black, green and orange circles show the magnetic length for three 
points in the membrane marked by a dot in the same colour. 
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